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Abstract: The modified evolution equation for parton distributions of Dokshitzer, March- 
esini and Salam is extended to non-singlet Deep Inelastic Scattering coefficient functions 
and the physical evolution kernels which govern their scaling violation. Considering the 
x — > 1 limit, it is found that the leading next-to-eikonal logarithmic contributions to the 
physical kernels at any loop order can be expressed in term of the one-loop cusp anomalous 
dimension, a result which can presumably be extended to all orders in (1 — x), and has 
eluded so far threshold resummation. Similar results are shown to hold for fragmenta- 
tion functions in semi-inclusive e + e~ annihilation. Gribov-Lipatov relation is found to be 
satisfied by the leading logarithmic part of the modified physical evolution kernels. 
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1. Introduction 

There has been recently renewed interest [1-9] in threshold resummation of "next-to- 
eikonal" logarithmically enhanced terms which are suppressed by some power of the gluon 
energy (1 — x) for x — > 1 in momentum space (or by some power of 1/N, N — > oo 
in moment space). In particular, in [2,5-9] this question has been investigated at the 
level of "physical evolution kernels" which control the scaling violation of (non-singlet) 
structure functions. The scale-dependence of the Deep Inelastic Scattering (DIS) coef- 
ficient function C?,{x, Q 2 , n 2 F ) corresponding to the flavor non-singlet F2(x,Q 2 ) structure 
function^ (x, Q 2 )/x = Ciix, Q 2 , ^ 2 p)®<l2,ns{x, fAp), where q2,ns(x, fi F ) is the corresponding 
quark distribution) can be expressed in terms of C2(x,Q 2 , /xp) itself, yielding the following 
"physical" evolution equation (see e.g. Refs. [10-14]): 

(i.i) 

where fiF is the factorization scale (I assume for defmitness the MS factorization scheme is 
used). K(x, a s (Q 2 )) is the momentum space physical evolution kernel, or physical anoma- 
lous dimension; it is independent of the factorization scale and renormalization-scheme 
invariant. 
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In [15], the result for the leading contribution to this quantity in the x — > 1 limit was 
derived, which resums all logarithms at the leading eikonal level, and nicely summarizes 
analytically in momentum space the standard results [16, 17] of threshold resummation: 

K(x, a s (Q 2 )) ~ + Bf IS (a s (Q 2 )) 5(1 - x) , (1.2) 

where r = — - (with rQ 2 = W 2 the final state "jet" mass), Bjf* IS (a s ) is related to the the 
quark form factor, and J(Q 2 ), the "physical Sudakov anomalous dimension" (a renormal- 
ization scheme invariant quantity), is given by: 



(1.3) 

oo 

1=1 

In eq.(1.3), 

oo 

A{a 8 ) = ^A i a i 8 
i=i 

is the universal "cusp" anomalous dimension [18] (see also [19]), with a s 
coupling, 

P{a s ) = = "A) a 2 -fha* s - ft *\ + ... (1.5) 

is the beta function (with (3q = — \ n f) an( i 

oo 

B(a s ) = Y J B l a\ (1.6) 

i=\ 

is the usual final state "jet function" anomalous dimension. It should be noted that j\ = A\ 
(the one loop cusp anomalous dimension), and also that both A(a s ) and B(a s ) (in contrast 
to J{Q 2 )) are renormalization scheme-dependent quantities. The renormalization group 
invariance of J(Q 2 ) yields the standard relation: 



J(Q 2 ) = A(a s (Q 2 )) + 



dB (a s (Q 2 )) 



dlnQ 2 



A(a s {Q 2 ))+p(a s {Q 2 )) 



dB(a s (Q 2 )) _ 



da q 



(1.4) 



= % the MS 



J ((1 - x)Q 2 ) = ji a s + a 2 [-j 1 p L x + j 2 } (1.7) 
+ a 3 s [j^ 2 L 2 x - fafa + 2j 2 fo)L x + j 3 ] 

+ a A s [-hPlLl + (|ji/3iA) + ZhPl)L 2 x - (jifo + 2j 2 /3i + 3j 3 fo)L x + j 4 ] + ... , 

where L x = ln(l— x) and a s = a s (Q 2 ), from which the structure of all the eikonal logarithms 

in K(x,a s (Q 1 )) can be derived. A term like arising from — ^— — - in eq.(1.2) must be 
interpreted as usual as a standard +-distribution. the eikonal logarithms are thus 
absorbed into the single scale (1 — x)Q 2 (see also [20-22] and section VI-E in [11]). 
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However, no analogous result holds [6] at the next-to-eikonal level (except [2] at large- 
fio). In this note, I show that the leading next-to-eikonal logarithmic contributions to the 
physical evolution kernel at a given order in a s can actually be determined in term of lower 
order leading eikonal coefficients, representing the first step towards threshold resummation 
at the next-to-eikonal level. This result is obtained by extending the approach of [23, 24] 
(which deals with parton distributions) to the DIS coefficient functions themselves. 



2. The modified physical kernel 

I consider the class of modified physical evolution equations: 



dC 2 {x,Q 2 ^ 2 F ) = /•' dz h ii)2 iw> / ,^ n2/ .\ 

J X 



■K(z,a s (Q 2 ),X)C 2 (x/z,Q 2 /z\^ 2 F ) , (2.1) 



dlnQ 2 J x z 

where for book-keeping purposes I introduced the parameter A, which shall eventually be 
set to its physically meaningful value A = 1, in straightforward analogy to the modified 
evolution equation for parton distributions of [24]. I note that K(x,a s ,X = 0) = K(x,a s ), 
the 'standard' physical evolution kernel. Eq.(2.1) allows to determine K(x,a s , X) given 
K(x,a s ) (or vice-versa). Indeed, expanding C2{y, Q 2 /z x , fx 2 F ) around z = 1, keeping the 
other two variables fixed, and reporting into eq.(2.1), one easily derives the following rela- 
tion between K(x,a s ,X) and K(x,a s ): 



K(x,a s ) = K(x, a s , A) — A [Inx K(x, a s , A)] K(x, a s ) (2-2) 
+y[ln 2 x K{x,a s , A)] <g> \P(a 8 ) dK ^ a ** + K{x, a s ) ® K(x, a s )} + ... , 

where only terms with a single overall factor of A need actually to be kept up to next-to- 
eikonal order, since one can check terms with more factors of A, which are associated to 
more factors of Inx, are not relevant to determine the next-to-eikonal logarithms in the 
physical kernel. In the rest of the paper (except section 4) I shall therefore simply use: 

K(x, a s ) = K(x, a s , A) — A [In x K(x, a s , A)] <g> K(x, a s ) + ... . (2-3) 
Eq.(2.3) can be solved perturbatively. Setting: 

K(x, a 8 , A) = K (x, X)a s + K^x, X)a 2 + K 2 (x, X)a 3 s + K 3 (x, X)a 3 s + ... (2.4) 
(and similarly for K(x,a s )), one gets: 



K (x, A) = K (x) (2.5) 
K^x, A) = K^x) + X[\nx K (x)} ® K (x) 

K 2 (x, A) = K 2 (x) + X{[lnx K^x)] ® K (x) + [Inx K (x)] ® K^x)} + ... 
K 3 (x, A) = K 3 (x) + A{[lna; K 2 (x)} ® K (x) + [Inx K^x)} ® K^x) 
+ [\nx K (x)] <g) K 2 (x)} + ... . 
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The Ki(x) , s are determined in term of splitting functions and coefficient functions as follows 
[14]: 



K (x) = Po(x) (2.6) 

ifi(x) = Pi(x)-A) ci(x) 

^ 2 (x) = P 2 (x) - /3i ci(x) - /3 (2c 2 (x) - cf (x)) 

tf 3 (s) = P 3 (x) " & ci(x) - /3i(2c 2 (x) - cf 2 (x)) - /3 (3c 3 (x) - 3c 2 (x) ® ci(x) + cf 3 (x)) , 

where Pj(x) are the standard (i + l)-loop splitting functions, q(x) are the i-loop coefficient 
functions, and cf 2 (x) = ci(x) <g> ci(x), etc.... 

Consider now the x — > 1 limit. The one- loop splitting function is given by [25]: 

P (x) = A lPqq {x) + B{ 6(1 - x) , (2.7) 

with Ai = 4Cf, and 1 : 

= - 1 + ^ - *) = + \ {l - x) = \ + \ {l - x) ■ (2 - 8) 

Moreover, at the next-to-eikonal level we have, dropping from now on 5 function contribu- 
tions: 



with [26]: 



P 1 (x) = ^ + C 2 L X + D 2 + ... , (2.9) 



C 2 = A\. (2.10) 



Also: 

Cl{x) = CllL ^ +Cl ° + bll L x + b 10 + ... (2.11) 

with c\\ = Ai = 4Cp, bn = 0. From eq.(2.6) one can derive [6,7] the following expansions 
for x — > 1: 



K (x) = P (x) = — + h w + ... (2.12) 

r 

, x ^lAr + &20 . , r . , 
K\[X) = h fol-Lx + /l20 + ••• 

r 

k 32 Ll + A^iL^ + fe 30 



K 2 {x) = -2—2 ^ 2^ + h 32 L^ + + ^ + ... 

r 

fc 43 £3 + ^42^ 2 + hiL x + fc 40 3 2 

if 3 (a;) = h /i43^ + h i2 L x + /i^Lj. + h m + .. 



1 Pqq(a;) is defined to be 1/2 the corresponding function in [7]. 



- 4 - 



3. Leading next-to-eikonal logarithms 

3.1 Two loop kernel 

From eq.(2.6), (2.7), (2.9) and (2.11) one deduces: k 10 = A ± , h w = 0, and k 21 = -fioA 1 , 
h 2 \ = C 2 . Then eq.(2.5) yields for x ->■ 1: 

K (x,X) = P (x) (3.1) 
Ki{Xt A) = k 21 L x + k 20 + _ A ^ + o{L x) 

Now 

MA) = h 21 - A 4, = C 2 - A A? = (1 - \)A\ . (3.2) 

Thus, setting A = 1, one finds that the leading next-to-eikonal logarithm in Ki(x, A = 1) 
vanishes, yielding the relation: 

h 21 = k\ = 16C£ , (3.3) 

which is correct [6,7]. This finding is not surprising: up to two loop, the leading next-to- 
eikonal logarithm is contributed only by the splitting function, since b\\ = (e.g. h 2 \ = C 2 ), 
and one effectively recovers the result (eq.(2.10)) holding [24] for the two loop splitting 
function. The situation however changes drastically at three loop, where the leading next- 
to-eikonal logarithm is contributed by the coefficient function rather then the splitting 
function, and the crucial question is whether the leading next-to-eikonal logarithm still 
vanishes for A = 1. 

3.2 Three loop kernel 

Eq.(2.5) yields for x 1: 

K 2 (x, A) = h2Ll + hlL * + ho + (h 32 - \h 21 k w )Ll + 0(L X ) . (3.4) 
Requiring /132(A), the coefficient of the O(L^) term, to vanish for A = 1 predicts: 

= |feifeio = ~\^A\ = -24/3 C| , (3.5) 

which is indeed the correct [6, 7] value. I stress that this result is not a consequence of the 
relation [24,27,28] C 3 = 2A 1 A 2 for P 2 (x). Indeed it is well-known [29] that the P,(a;)'s, 
and in particular P 2 (x), have only a single next-to-eikonal logarithm: 

P 2 (x) = — + C 3 L X + D 3 + ... , (3.6) 
r 

and thus P 2 (x) cannot contribute to the double logarithm in K 2 {x). Rather, h 32 is con- 
tributed by the coefficient functions in eq.(2.6), and eq.(3.5) yields a prediction for the 
0(L X ) term in c 2 (x). 
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3.3 Four loop kernel 

Eq.(2.5) yields for x — Y 1: 



fc 43 L3 + /c 42 L2 + k 41 L x + /c 40 
r 

+[h 43 - \{h w k 32 + \k 2 21 )]Ll + 0{L 2 X ) , (3.7) 

where k 32 = ^4i/?o (consistently with eq.(1.7)). Requiring /143(A), the coefficient of the 
0(L X ) term, to vanish for A = 1 predicts: 

&43 = ^10*32 + ^Ii = y$A? = y/3oCf , (3-8) 
which is again the correct [6, 7] value. 

3.4 Five loop kernel 

One can similarly predict the leading next-to-eikonal logarithm in the five loop physical 
kernel (which depends on the four loop coefficient function). Using eq.(2.3), the coefficient 
of the 0(L*) term in K 4 (x, A) is found to be given by: 

5 5 

/154(A) = /154 - A(-/ci &43 + gfc 2 i/c 32 ) , (3.9) 

where £43 = — A\/3q (again consistent with eq.(1.7)). Requiring this coefficient to vanish 
for A = 1 predicts 2 : 

5 5 25 100 

h 54 = -k w k 43 + -k 2l k 32 = -r^PlAl = -—ffiCl . (3.10) 

3.5 All-order relations 

Defining moments by 

K(N,a s )= [ dxx N - 1 K(x,a s ) , (3.11) 
Jo 

eq.(2.3) yields in moment space: 

K(N,g s )= K ( N > a »V , (3.12) 
V ' l + XK(N,a s ,X) K ' 

where / = df/dN. Assuming the leading next-to-eikonal logarithms vanish to all orders 
(in an expansion in 1/r) for A = 1, i.e. that /ij+i,i(A = 1) = for i > 0, one can derive [30] 
from eq.(3.12) the resummation formula: 

00 . . 

J> +M LX +1 = Hl + a s f3 L x ) , (3.13) 

^— ' po 1 + a s p L x 

1=0 

which correctly reproduces the results in the previous subsections. 



2 Hence £dis 4 = in the notation of [7]. 
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One can further show [30] that the moment space functional relation which accounts for 
leading logarithms at all orders in (1 — x) is: 



K(N, a s ) = K[N - XK(N, a s ), a s ,X] . (3.14) 

Eq.(3.12) results from expanding the right hand side of eq.(3.14) to first order in AiV = 
XK(N,a s ). It is interesting that eq.(3.14) is identical to the functional relation 3 obtained 
[28, 32] for the splitting functions in the conformal limit (where the splitting functions 
coincide with the -fQ's). 



4. Leading next-to-next-to-eikonal logarithms 

It can be checked [30] that similar methods allow to predict using eq.(2.2) the leading 
logarithmic contributions at the next-to-next-to-eikonal level, i.e. the coefficient of the 
(1 — x)U x term in Ki(x). The crucial new point, however, is that the leading term in the 
eikonal expansion has to be defined in term of the one-loop splitting function prefactor 
p qq {x) (eq.(2.8)), instead of 1/r as in eq.(2.12). Namely, keeping only leading logarithms 
at each eikonal order, the predicted coefficients (j = i + 1, i > 0) are defined 4 by: 



Ki(x)\ LL = Ll\p qq (x) kji + hji + (1 - x)f 3l + (1 - xfg 3l + - x) 3 )] . (4.1) 
Eq.(2.2) yields the corresponding fji(X) coefficients in Ki(x,X): 

/ 2 c i(A) = Zli + A^o (4.2) 



/43(A) = /I3 - A^ - ^k w k 32 + -(/121 - ^Ofoi + ^10^32) + A 2 /ci fc 2 i , 

where I used that hi = 0, and one should note the presence of contributions quadratic in 
A. Assuming the /Jj(A)'s vanish for A = 1, one thus derives the relations (with /f = 0): 

fli = ~\k 2 w = -8C| (4.3) 
/ 3 C 2 = ~\k W k 2 i + k w h 21 - ^kf = + 32C 3 F 

/Is = ~\k w ks 2 + ^{h 21 - h 21 )k 21 + k 10 h 32 - k 2 l0 k 21 = ~y c f^ ~ 64C£/3 , 

which are seen to be correct using eq.(3.26) in [7]. The latter equation also makes it likely 
that similar leading logarithmic predictions can be obtained to any order in (1 — x), using 



3 A similar functional relation has been obtained in a different context in [31]. 

4 The motivation for the superscript "c" (for "classical") shall be clarified in the Conclusion section. 



-7- 



the same prefactor p qq {x) as in eq.(4.1) to define the leading term in the eikonal expansion. 
Indeed, one derives for instance [30] the 0((l — x)) 2 coefficients in eq.(4.1) (with gi$ = 0): 



921 = \k\o = fc* (4.4) 
532 = gfciofoi + (/21 + 2^21 + ^k 2 i)ho = 2^10^21 = -8Cp(3 , 

which are correct [7]. I note that /fi and 521 coincide (like /i2i) with the splitting functions 
contributions. 



5. Fragmentation functions in e + e annihilation 

Similar results hold for physical evolution kernels associated to fragmentation functions 
in semi-inclusive e + e~ annihilation (SIA), provided one sets A = — 1 in the analogue of 
eq.(2.1): 

^tinQ^ = j x 1 ^KT(z,a s (Q%\)C T (x/z,Q 2 /z\ti) , (5.1) 

where Ct denotes a generic non-singlet SIA coefficient function. I first note that threshold 
resummation in this case [33] leads at the leading eikonal level to an equation similar to 
eq.(1.2): 

K T (x,a s (Q 2 )) ~ J{{ \~_ X) x Q2) +Bl IA {a s {Q 2 ))5{l-x) , (5.2) 

where x should now be identified to Feynman-x rather then Bjorken-x, and I used the 
results of [34] which imply that the "physical Sudakov anomalous dimension" J(Q 2 ) is the 
same for structure and fragmentation functions. The statement above eq.(5.1) then follows 
from the following two observations: 

i) The predictions in eq.(3.3), (3.5), (3.8) and (3.10) depend only upon coefficients of leading 
eikonal logarithms in the physical evolution kernels. 

ii) Eq.(3.26) in [7] shows that the latter coefficients are identical for deep-inelastic struc- 
ture functions and for e + e _ fragmentation functions (consistently with the remark below 
eq.(5.2)), but that the coefficients of the leading next-to- eikonal logarithms are equal only 
up to a sign change (in an expansion in 1/r) between deep- inelastic structure functions 
and fragmentation functions. 

One deduces the resummation formula (j = i + 1): 

00 A A 

~^ Po 1 + a s p ^x 
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6. Conclusion 

A modified 5 evolution equation for DIS non-singlet structure functions, analoguous to 
the one used in [24] for parton distributions, but which deals with the physical scaling 
violation and coefficient functions, has been proposed. It allows to relate the leading next- 
to-eikonal logarithmic contributions in the momentum space physical evolution kernel to 
coefficients of leading eikonal logarithms at lower loop order (depending only upon the 
one- loop cusp anomalous dimension A\), which represents the first step towards threshold 
resummation at the next-to-eikonal level. This result also explains the observed [6, 7] 
universality of the leading next-to-eikonal logarithmic contributions to the physical kernels 
of the various non-singlet structure functions, linking them to the known [35] universality 
of the eikonal contributions. Similar results hold at the next-to-next-to-eikonal level with a 
proper definition of the leading eikonal piece, and can presumably be extended to leading 
logarithmic contributions at all orders in (1 — x). Analogous results are obtained for 
fragmentation functions in semi-inclusive e + e~ annihilation. 

One may ask to what extent the success of the present approach may be attributed, 
as suggested in [24, 32] for the splitting functions case, to the classical nature [36] of soft 
radiation. In fact, the main result of this paper for the (modified) DIS physical evolution 
kernel can be summarized (barring the 5-function contribution) by the following equation: 



K(x,a s , A = 1) ~ 



X 1 /-. X 



J ((1 — x)Q 2 ) + subleading logarithms , (6.1) 



where the second term (the "subleading logarithms") is contributed by all powers in (1 — x) 
except the leading eikonal one. The first term in eq.(6.1) accounts for the leading logarith- 
mic contributions to the modified kernel (together with some subleading logarithms) to all 
powers in (1 — x) at any given loop order, and implies leading logarithmic contributions 
are actually absent beyond 0(1 — x) power. This term has the remarkable effective one- 
loop splitting function form 4Cf a p hys ((1 — X )Q 2 )) Pqq{ x )-, with the "physical coupling" 

As pointed out in [32], the part of the one-loop prefactor (eq.(2.8)) should be 
interpreted as corresponding to universal classical radiation, a QCD manifestation of the 
Low-Burnett-Kroll theorem [36], while the 1 — x part represents a genuine quantum con- 
tribution. Now, it is clear that at the next-to-eikonal level, the 1 — x part of the prefactor 
is irrelevant: only the "classical" 1/r part is required to separate those leading logarithms 
in the standard (A = 0) physical evolution kernel which are correctly predicted in the 
present approach (the hji in eq.(4.1)), hence "inherited" in the sense of [32], from the "pri- 
mordial" ones (those which at each loop order carry the same color factors as the leading 
0(1/(1 — x)) eikonal logarithms, and can thus be absorbed into the definition of the leading 
term). However, it appears from the results of section 4 that, at next-to-next-to-eikonal 
level, the full one-loop prefactor has to be used into the definition of the leading term 



5 Evolution equations involving similar kinematical rescaling factors have been suggested in the past: see 
e.g. eq.(3.2) in [20]. 
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to properly isolate the "inherited" next-to-next-to-eikonal logarithms (the in eq.(4.1)). 
Moreover, although the "inherited" are purely "classical" (like the hji), the "inherited" 
gji at the 0((1 — x) 2 ) level are a mixture of "quantum" and "classical". Indeed, setting 
/io = 5^10 and /21 = 5^21 (the "quantum parts" of the 0(1 — x) coefficients), one finds 

521 = 521 + 521 > with 521 = hofto = l k 10 and 521 = -^10; and 532 = 532 + 532> with 

532 = ^21/10 + ^10/21 = jhiho and g c m = -jfoi&io- In both cases c/£ = -\g%, which 
shows the "inherited" gji coefficients are actually dominantly "quantum". 

It can be further checked [30] that the very same first term in eq.(6.1) also accounts for 
the leading logarithmic contributions to the A = — 1 modified SIA physical evolution kernel 
to all powers in (1 — x), which implies that the leading logarithmic parts of the modified DIS 
and SIA physical evolution kernels satisfy Gribov-Lipatov relation [37], namely we have: 

K(x, a s , A = 1)| LL = K T (x, a 8 , A = -1)| LL = p qq (x) J ((1 - x)Q 2 ) | LL , (6.2) 

where J ((1 - x)Q 2 ) | LL = A (a s ((l - x)Q 2 )) | LL = 1+ ^q ( 2)2l x is the leadm S logarithmic 
contribution to eq.(1.7). Indeed, once transformed back to the standard (A = 0) physical 
kernels, eq.(6.2) is consistent with eq.(3.26) in [7] at least to next-to-next-to-eikonal order, 
and is probably correct to all orders in (1 — x) (with identically vanishing contributions 
beyond 0(1 — x) order). On the other hand, contrary to the splitting functions case where 
it has been checked up to three loops [32,38], a full Gribov-Lipatov relation K(x, a s ,\ = 
1) = Kt{x, a s , X = — 1) does not seem to hold for subleading logarithms beyond the leading 
cikonal level. 

The resummation of the subleading logarithmic contributions at next-to-eikonal order 
in eq.(6.1), not adressed here, remains an open issue: the present method does not work 
for them, except in the conformal limit, where one recovers the results of [24]. 

Acknowledgements I thank G. Marchesini for stimulating discussions and hospitality at 
the University of Milano-Bicocca where part of this paper has been written. I also wish to 
thank the referees for constructive suggestions. 
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